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At the present time there is a lively debate within the nuclear community concerning the relevance 
of quark degrees of freedom in understanding nuclear structure. We outline the key issues and review 
the impressive progress made recently within the framework of the quark-meson coupling model. 
In particular, we explain in quite general terms how the modification of the internal structure 
of hadrons in-medium leads naturally to three- and four-body forces, or equivalently, to density 
dependent effective interactions. 
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' I. INTRODUCTION 

<n : 

We have a fundamental theory of the strong interaction, namely Quantum Chromodynamics (QCD), the math- 

1 ^ • ematical beauty of which, in combination with its phenomenological successes, has convinced most physicists that 

it must be correct. Its predictions have been accurately confirmed in the short-distance (or high energy) regime 
of "asymptotic freedom". In the opposite limit of long distances (corresponding to quark confinement), it will be 
extensively tested over the next decade through advances in lattice QCD combined with new experimental facilities, 
H ■ such as the 12 GeV upgrade at Jefferson Lab and the new hadronic capabilities at J-PARC and FAIR. One of the 
^ ' most compelling challenges in modern nuclear physics is to relate the properties of nuclear matter, from finite nuclei 
C^ - , to neutron stars, to the underlying quark and gluon structure of matter and QCD itself. 

Over the past 50 years the standard approach to understanding nuclear structure has involved the application of non- 
relativistic many-body theory, based upon nucleon-nucleon (NN) potentials fit to experimental two-nucleon scattering 
data. Amongst the potentials used at various stages we mention the Paris potential jl[, with its intermediate range 
structure determined through dispersion relations and hence almost model independent. Other approaches which are 
still widely used are the one-boson-exchange forces, such as the various Bonn [2[ and Nijmegen 0] potentials. In its 
' modern form this approach is best represented by the Argonne 18 potential [1], supplemented by a phenomenological 
three-body force Q. This has been used in combination with Green function monte-carlo methods to calculate the 
1 1 ■ energy levels of light nuclei. 

These models of the nuclear forces have no direct connection to QCD, in that they operate at the purely hadronic 
level. Another recent approach, which preserves the chiral symmetry of QCD Q, is based on effective field theory 
(EFT). Such calculations are just now achieving the phenomenological success in fitting NN data that is characteristic 
of the Bonn and Nijmegen potentials mentioned above - with a similar number of fitting parameters, typically between 
20 and 30. The apparatus of EFT is also being exploited to generate the corresponding three-body forces Q, so that 
one can tackle the energy levels of light nuclei. While the use of EFT permits one to preserve the chiral symmetry of 
QCD, this is a pale reflection of the full power of QCD itself. Indeed, given the extensive knowledge of chiral symmetry 
at the time, one could have carried out the EFT program for the NN force in the 1960s, before the discovery of QCD. 

Neither the usual treatments of nuclear structure outlined above, nor the usual EFT approach, have yet incorporated 
the requirements of special relativity. We mentioned earlier the results of the Paris group's dispersion relation analysis 
of the nature of the intermediate range NN force. The intermediate range NN force is unambiguously dominated by 
two-pion exchange with a Lorentz scalar and isoscalar character. Quantum Hadrodynamics (QHD) exploited this by 
formulating a description of nuclear matter in which Lorentz scalar attraction (represented by a a meson) competed 
with Lorentz vector repulsion (represented by an ui meson) Q. At the saturation density of nuclear matter, pQ, the 
mean scalar field strength in QHD was of order 400 MeV. A similar result was found in a relativistic Brueckner 
Hartree-Fock calculation using a boson exchange potential @ . Given such a huge scalar field, with a strength almost 
half of the mass of the nucleon itself, one must anticipate that it will have a significant effect on the internal structure 
of the nucleons making up the nuclear matter. In particular, we recall that the typical energy associated with the 
excitation of internal degrees of freedom in the nucleon is only 300 MeV (for spin excitations) to 500 MeV (for orbital 
excitations). 

Let us briefly summarize. Any serious treatment of the model independent Lorentz scalar nature of the intermediate 
range NN force leads to the conclusion that the bound nucleons experience a scalar potential of a strength up to one 
half of the mass of the nucleon itself - a strength comparable with the typical internal excitation energy of the 
nucleon. This immediately suggests that the internal structure of the nucleon should indeed play a crucial role in 
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nuclear structure. The popular prejudice which leads many to neglect this role is generally based on the observation 
that the nucleon binding energy is very small in comparison with its mass. However this small binding actually 
arises from a cancellation between the scalar and the vector potential in-medium, while the structure modification is 
expected to arise primarily through the scalar potential. 

In order to anticipate the possible consequences of this very general observation, we turn to atomic and molecular 
systems where we have text book experience. We know that when an atom is subjected to a strong electric field, its 
electron structure will rearrange in order to oppose the applied field. This change in the internal structure of the 
atom, at least if one is concerned solely with describing the energy of the system, can be described in terms of an 
electric polarizability. In particular, the energy of the system has a term quadratic in the applied electric field, with 
the coefficient being (one half of) the electric polarizability. Exactly the same thing happens if we apply a magnetic 
field, with the coefficient of the term quadratic in the applied magnetic field involving the magnetic polarizability. 

Turning to the nucleon itself, we know that applied electric and magnetic fields alter its internal structure, giving 
rise to electric and magnetic polarizabilities and even, in sophisticated electron scattering experiments, the so-called 
generalized polarizabilities [13, Given this background it is remarkable that it has taken so long for nuclear 

physicists to pay serious attention to the response of the nucleon to an applied scalar field, especially its scalar 
polarizability. 

Taking the experience with atomic physics as our lead, we would naturally expect that a nucleon embedded in 
matter would have an energy with a non-linear dependence on the mean scalar field: 

M* N = M N - g a a + ^(g a a) 2 (1) 

where d is the scalar polarizability. In QHD such a term is, of course, absent - although later phenomenological 
developments involving a self-interaction of the scalar field can be rewritten (using a redefinition of the scalar field) 
in just this way. Viewed in this way the scalar polarizability can be seen as giving a very natural explanation of such 
a self-interaction. In the first examination of nuclear matter from this point of view, using the MIT bag model to 
describe the quark structure of the bound nucleons, Guichon found exactly such a behavior, with d = +0.22R and 
R the bag radius [ll|. Since that time, this model, which is known as the quark-meson coupling (QMC) model, has 
been extended to describe finite nuclei [Hj]. It has also been widely applied to the extremely interesting problem of 
how the structure of hadrons in general might be expected to change in-medium [l4l flHl - applications which have 
included the possibility of rj- or cj-nucleus bound states (lal. the structure of hypernuclei |l7| and the modification of 
the electric and magnetic form factors of bound protons [l8j] , as well as their structure functions . 

Rather than reviewing these many applications, which have recently been the subject of a major review we 
return to the atomic and molecular analogy for further insights which, as we shall see, are crucial to understanding 
the structure of atomic nuclei. When two atoms approach each other they mutually modify their internal structure. 
This implies that when a third atom approaches it does not see a pair of bare atoms. Consequently, the total energy 
of the system will not be the sum of the individual pair-wise interactions (potentials) - that is: 

Vtot + V 12 + V 23 + Vi 3 (2) 

and the difference is, by definition, a "three-body force" , V123. By analogy, we anticipate that the scalar polarizability 
of the nucleon will naturally lead to many-body forces in-medium. 

In the next section we briefly review the exploration of the origin of many-body forces (or equivalently density 
dependent effective interactions), arising as a consequence of the modification of the internal structure of the nucleon 
in- medium, within the framework of QMC. We shall see that it is indeed possible, starting from the quark level, to 
generate realistic effective interactions of the Skyrme type ■ 



II. EFFECTIVE FORCES OF THE SKYRME TYPE DERIVED FROM THE QUARK LEVEL 

In 2004 Guichon and Thomas [2l| showed that by expanding about < a >= one could derive an effective force 
of the Skyrme type (widely used in nuclear structure calculations) in which the local two-body effective interaction 
was supplemented by three- and four-body forces, proportional to d and d 2 , respectively. A comparison between the 
various coefficients in the Skill force with those derived from the QMC model showed agreement at typically the 10% 
level. This remarkable result provided the first, direct microscopic connection between a commonly used effective 
nuclear force and the underlying degrees of freedom of QCD. 

As important as the results of Ref. [2l[ were, many important issues were left unaddressed: 

• The relatively large scalar field found in nuclear matter means that one would really like to remove the need to 
expand about a = 0. 
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• Modern Skyrme forces tend to include density dependent coefficients (fit to appropriate nuclear data), rather 
than being density independent with an explicit three-body term. Hence one would really like to rewrite the 
energy functional for QMC in terms of a density dependent effective force. 

• A crucial issue at higher densities, when for example one is calculating the equation of state (EoS) for application 
to neutron stars, concerns the density at which hyperons appear in matter in /3-equilibrium. For this purpose 
one would very much like to derive effective, density dependent effective hyperon-nucleon and hyperon-hyperon 
forces. Because one starts from the quark level, QMC is ideally suited for this purpose. 

All of these issues were successfully resolved by the recent formulation of Guichon et al. [12] • The resulting energy 
functional in its full relativistic form was too difficult to use in finite nuclei. However, in that case a non relativistic 
expansion was allowed and that provided an effective force with a novel density dependence - the coefficients being 
rational functions of the local density. For example, the central pieces of the corresponding Hamiltonian 



<H(r) >=pM+^+?A)+H 



took the form: 
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This appears rather different from commonly used density dependent Skyrme forces such as SkM, which includes a 
dependence on a fractional power of the density. For comparison, we show for SkM the central pieces corresponding 
to those that were given for QMC in Eq. (g]): 



P U 3 (2 p 2 - p n 2 - p p 2 ) t (p 2 (2 + X ) - (1+2X0) (Pn 2 +Pp 2 )) 
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In order to make a quantitative comparison between these two forms of effective potential the SkM form was fit to 
the QMC form over the density range (0,po)- The comparison between the usual phenomenological values of the 
constants appearing in the SkM force and those fitted to QMC showed a level of agreement at the 10-20% level. This 
is reasonable but by no means impressive. On the other hand, since the functional forms are quite different it seemed 
appropriate to compare the predictions of the two effective forces for real nuclear properties. Accordingly the effective 
force derived from the QMC model was employed within the Hartree-Fock framework to calculate the properties of 
closed shell, finite nuclei. The result was extremely gratifying, with the level of agreement between theory and data 
(where available and phenomenological Skyrme-type forces where data was unavailable) being very good. In Fig. [1] we 
show a comparison between the charge densities of selected closed shell nuclei calculated within QMC in comparison 
with experimental data, as well as the widely used Sly4 force [23| . Clearly the level of agreement between the two 
theories is remarkably good. Given that the QMC calculation had no parameters adjusted to reproduce the properties 
of finite nuclei, the agreement with the experimental charge densities is also very satisfactory. 

As there is no model independent experimental determination of the neutron densities we show, in Fig. [51 a 
comparison between the neutron densities produced by Sly4 and by the QMC model. Once again, it is reassuring in 
terms of assessing the quality of the effective interaction derived from the QMC model that they are so close. 

Another key test of the model is its ability to describe the dependence of the spin-orbit splittings on both mass 
number and isospin. At first glance this looks like an area in which QMC might fail, in view of the strong isospin 
dependence of the effective force derived from QMC. However, as we see in Table Q] (taken from Ref. [13]), this worry 
is in fact misplaced. Once the interaction is used, within the Hartree-Fock framework, to self-consistently determine 
the nuclear structure there is agreement between the predictions of QMC and experiment across the full range of 
nuclear mass number for both protons and neutrons. 



III. CONCLUDING REMARKS 



In the limited space available we have been able to introduce just the main ideas of these modern QMC based 
calculations. The results obtained with the density dependent effective force are particularly impressive. However, 
they are based upon the non-relativistic approximation, which cannot be expected to be reliable much beyond po- 
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FIG. 1: Proton densities calculated using the density dependent effective interaction derived from the QMC model compared 
with experiment and with the predictions of the Skyrme Sly4 force - from Ref. [%| . 
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FIG. 2: Neutron densities calculated using the density dependent effective interaction derived from the QMC model compared 
with the predictions of the Skyrme Sly4 force - from Ref. [22l [. 
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Neutrons (exp) 


Neutrons (th) 


Protons (exp) 


Protons (th) 


W 0, lpi/2 ~ lp 3 /2 


6.10 


6.01 


6.3 




40 Ca, lds/a - ld 5/2 


6.15 


6.41 


6.00 


6.24 


48 Ca, ld 3/2 - lds/a 


6.05 (Sly4) 


5.64 


6.06 (Sly4) 


5.59 


208 Pb, 2d 3/2 - 2d 5/2 


2.15 (Sly4) 


2.04 


1.87 (Sly4) 


1.74 



TABLE I: Values of the spin-orbit splitting for selected nuclear levels calculated from the QMC model, in comparison with the 
corresponding experimental values, where known. As they are not so well known in the case of 48 Ca and 20S Pb, there we give 
the values corresponding to the Skyrme Sly4 force. 
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FIG. 3: QMC predictions for the pressure as a function of density in comparison with the allowed region for isoscalar matter 
deduced from the constraints of heavy- ion flow data [25| . 

Indeed, as shown by Guichon et al. [22], it leads to an error of almost 50% in the velocity of sound in the nuclear 
medium at just 2pQ. Thus, in studying dense matter, one must return to using the original relativistic Hartree- 
Fock formulation. This is feasible in uniform matter [2J]. As an illustration of our results, in Figure [3] we show 
our predictions for the pressure of symmetric nuclear matter as a function of baryonic density. The labels QMC700 
and QMC7r4 refer to two extreme versions of the model which give an incompressibility of 340 MeV and 256 MeV, 
respectively. Also shown in the figure is the experimental constraint deduced from the flow data in high energy nuclear 
collisions [25j ]. 

Another challenge faced by nuclear theory, as the density increases, is the possible entry of hyperons into the 
nuclear stew. Hyperon-nucleon forces are unfortunately poorly constrained by data, even in the case of A N, while 
in the case of E N and S N the situation is worse. Typical phenomenological forces rely on SU(3) symmetry, a hard 
path to tread when cancellations between various components of the force are so important. When it comes to using 
hypernuclear data to constrain the parameters of Skyrme-type forces we note that only a single E-hypernucleus has 
been confirmed and no S-hypernuclei. For A-hypernuclei the situation is a little better. It is little wonder that there 
is no consensus on the threshold density at which hyperons appear, nor even which one will appear first. As the 
appearance of hyperons will soften the nuclear equation of state and, for example, raise the critical density for the 
transition to quark matter, this is a vital issue. 

Because the QMC model starts at the quark level, it can can also be used to calculate density dependent, effective 
hyperon-nucleon or even hyperon-hyperon forces - with no additional parameters. The only hypothesis needed is 
that the nuclear fields do not couple to the strange quark. This is certainly a good approximation if the meson fields 
simulate multi-pion exchanges. It is also supported by the absence of spin-orbit splitting in A hypernuclei. As shown 
in Figure [4] a striking prediction of the model is that the thresholds for the appearance of strange particles do not 
follow the usual pattern where the E _ and A appear first. Because of the scalar polarizability of the baryons and 
the use of the Hartree-Fock approximation [24| (rather than the mean field approximation), the 3 and A appear first 
and at a rather high density (~ 3po). The E does not show up in the density domain pertinent to neutron stars. 
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n B [fm" 3 ] 

FIG. 4: The f3 equilibrium population of particles calculated for the models QMC700 and QMC7r4 - from Ref. [24| . 

The resulting equation of state is relatively stiff and produces a maximum neutron star mass of about 2.1-Mq, with a 
rather small central density (~ 6/9o). 

Finally, we note that the QMC model is built upon the MIT bag model. This is a fairly ancient and unsophisticated 
quark model and although all indications are that the details of the model are not so important one would like to 
do better. In particular, one would ideally like to have a quark model describing the structure of the hadrons which 
is covariant, in order to address the problem of nuclear structure functions - the famous EMC effect [26|, [2?| - with 
some confidence. 

In parallel with the developments concerning QMC, Bentz, Thomas and collaborators have made considerable 
progress with just such a model. They took the NJL model, modified using proper time regularization in order to 
simulate confinement. The model is therefore covariant, confining and respects chiral symmetry. Exactly as in QMC, 
the hadron structure in-medium was self-consistently modified (by solving the corresponding Faddeev equations). 
Also as in the QMC model the scalar polarizability led to the saturation of nuclear matter and, incidentally, solved 
the long-standing problem of chiral collapse suffered by the NJL model 28] . In application to the nuclear EMC effect 
the model is able to reproduce the data for unpolarized structure functions (29l. [30|. However, even more important, 
it makes a prediction characteristic of models of the QMC type in which the mean scalar field enhances the lower 
Dirac components of the confined quark wave functions. (We note that a similar calculation within the chiral quark 
soliton model [3l[ produces a qualitatively similar result.) That is, it predicts a polarized EMC effect roughly twice 
as big as the unpolarized effect. This will provide a critical test of such models. 

A final success of this approach has been to use the same underlying quark model to describe matter made of 
"nucleons" (albeit with modified internal structure), as well as matter made of quarks [32l [33| . Indeed a major 
discovery of Lawley and collaborators [Hf was that, once the effect of the pion cloud of the nucleon [34| is taken 
into account, one can use the same underlying, quark level Hamiltonian to consistently describe the hadronic, quark 
matter and superconducting quark matter phases of dense matter. This is a major step forward in a field where one 
is usually forced to employ completely different models for the different phases. Such calculations are of great interest 
in view of the difficulties in simulating supernova explosions, as well as because of the very different mass-radius 
relations that one typically finds for stars containing quark matter. 
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